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Abstract — CHOKe is a simple and stateless active queue man- 
agement (AQM) scheme. Apart from low operational overhead, 
a highly attractive property of CHOKe is that it can protect 
responsive TCP flows from unresponsive UDP flows. Particularly, 
previous works have proven that CHOKe is able to bound both 
bandwidth share and buffer share of (a possible aggregate) 
UDP traffic (flow) on a link. However, these studies consider, 
and pertain only to, a steady state where the queue reaches 
equilibrium in the presence of many (long-lived) TCP flows and 
an unresponsive UDP flow of fixed arrival rate. If the steady 
state conditions are perturbed, particularly when UDP traffic rate 
changes over time, it is unclear whether the protection property 
of CHOKe still holds. Indeed, it can be examined, for example, 
that when UDP rate suddenly becomes (i.e., flow stops), the 
unresponsive flow may assume close to full utilization in sub- 
RTT scales, potentially starving out the TCP flows. To explain 
this apparent discrepancy, this paper investigates CHOKe queue 
properties in a transient regime, which is the time period of tran- 
sition between two steady states of the queue, initiated when the 
rate of the unresponsive flow changes. Explicit expressions that 
characterize flow throughputs in transient regimes are derived. 
These results provide additional understanding of CHOKe, and 
give some explanation on its intriguing behavior in the transient 
regime. 

Index Terms— Queue Management, RED, CHOKe, TCP, Flow 
Protection 

1. Introduction 

A. Overview on Flow Protection 

Broadly speaking, there are two distinct but complementary 
ways to enforce flow fairness and protection in the Internet. 
Following the end-to-end architectural design principle of the 
Internet H), the more classical way has been via conges- 
tion control algorithms 0. These algorithms are typically 
implemented in the transport protocols (e.g., TCP) of end 
hosts. To ensure global fairness, such schemes require all users 
to adopt them and respond to network congestion properly. 
However, this requirement can hardly be met for at least 
two reasons. Firstly, there is no performance incentive to end 
users. This is because users who lack the congestion control 
algorithms, intentionally or otherwise, may end up with a 
lion share of bandwidth. Secondly, in order to meet real-time 
requirements, many applications do not implement congestion 
control. Hence, to protect responsive (e.g., TCP) flows from 
unresponsive (e.g., UDP) ones, solely relying on the end-to- 
end schemes can be unfair or risky, and it is necessary to 
introduce some mechanisms in the network. This motivates 
the second approach to fairness and protection. 

* "Centre for Quantifiable Quality of Service in Communication 
Systems, Centre of Excellence" appointed by The Research Council 
of Norway, funded b y the Research Council, NTNU and UNINETT. 
http://www.q2s.ntnu.no 



The second approach is provided through router mecha- 
nisms. Such a router mechanism can be either (1) per flow 
fair queueing (PFFQ) scheme, e.g., Weighted Fair Queueing 
(WFQ) O, or (2) queue management (QM) scheme, e.g., 
Random Early Detection (PvED) |4), 0. PFFQ schemes share 
link bandwidth among flows in a fair manner. Typically, they 
isolate flows into separate logical or physical FIFO queues and 
maintain flow-level state information. By building firewalls 
around heavy users, flow isolation protects well-behaved flows 
and enables performance guarantees to such flows J6). Nev- 
ertheless, the maintenance of per flow state and the dynamic 
management of complex queue structure are widely believed 
to be problematicQ for high speed implementation. 

Queue management (QM) schemes specialize in buffer 
allocation and are normally significantly simpler in design, 
typically with a single queue shared by all flows. Among QM 
schemes, RED J4), (5) is probably the most widely known 
scheme. It maintains an exponentially moving average queue 
size that indicates the level of congestion in the router. A 
congested RED router drops incoming packets with a certain 
probability dependent on the queue size. Since the dropping 
probabilities are applied globally to all flows, both high rate 
and low rate flows can be punished in equal measures. In 
fact, based on the nature of Internet flows (e.g., flow sizes, 
underlying transport protocol) and resultant differences in their 
responsiveness to congestion, the same ambient drop rate of 
RED can be more detrimental and highly unfair to some 
flows. To deal with this unfairness, more complex variants 
of RED, e.g., Flow RED (FRED) Oj], RED with Propor- 
tional Differentiation (RED-PD) 1121 . have been proposed to 
apply differential per flow drop rates. However, such schemes 
typically need to maintain partial flow state to be able to 
discriminate drop rates among flows. 

With no flow isolation, fairness afforded by a QM scheme is 
generally approximate. For example, a recent simulation study 
consisting entirely of UDP flows shows that the bandwidth 
allocated by FRED to congested UDP flows can differ by 
several factors |9] see Fig. 6]. This begs the question: How 
much of the bandwidth, or generally the shared resources, 
can an unresponsive / rogue flow steal in QM schemes? 

B. CHOKe 

A highly novel QM scheme which, unlike FRED or RED- 
PD, does not require flow state to be maintained in the router, 

'Proponents of per flow architecture argue that the number of flows 
requiring scheduling at any moment is a lot less than the total number 
traversing the router. Hence, it is possible to realize scalable implementation 
of the per flow architecture (71 1101 . 
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is CHOKe |[T3l . CHOKe is proposed to protect rate-adaptive 
(responsive) flows from unresponsive ones. It uses the recent 
admissions, i.e. packets currently queued in the buffer, to 
penalize the high bandwidth flows. It can be implemented 
by a few tweaks of the RED algorithm. Specifically, "when a 
packet arrives at a congested router, CHOKe draws a packet 
at random from the FIFO buffer and compares it with the 
arriving packet. If they both belong to the same flow, then 
they are both dropped; else the randomly chosen packet is left 
intact and the arriving packet is admitted into the buffer with 
a probability (based on RED) that depends on the level of 
congestion M31I ." 

A promising property of CHOKe is that it provides analyti- 
cally proven protection of responsive flows from an unrespon- 
sive flow at the congested router, which provides an answer to 
the aforementioned question. Specifically, in the presence of 
many flows, the following interesting and peculiar steady-state 
properties of CHOKe have been derived lfl4l - irT6l : 

• Limits: An unresponsive UDP flow cannot exceed certain 
limits in buffer share and link bandwidth share ifTBl . 
ifTBI . The maximum UDP bandwidth share is (e + 1) = 
26.9% of link capacity, and the maximum buffer share is 
50%. 

• Asymptotic property: As the UDP rate increases without 
bound, its buffer share can asymptotically reach 50% 
and queueing delay can be reduced by half but its link 
utilization drops to zero. 

• Spatial distribution: The spatial packet distribution in the 
queue can be highly nonuniform Q4). The probability 
of finding a packet belonging to a high rate flow in the 
queue diminishes dramatically as we move towards the 
head of the queue. Correspondingly, the flow distribution 
in queue is skewed with most packets of high rate flows 
found closer to queue tail while packets of low rates are 
found closer to queue head. 

C. Contribution of This Paper 

To the best of our knowledge, previous analytical studies 
on CHOKe lfT4l - lfT6l are restricted to the steady state where 
the traffic rate of the UDP flow is assumed constant. However, 
this assumption is too restrictive, limiting more in-depth un- 
derstanding of CHOKe. This paper studies CHOKe behavior 
in the face of dynamically changing UDP rates and at the 
same time generalizes the steady-state properties proved in 
the earlier works. Particularly, this paper investigates CHOKe 
queue properties in a transient regime, which is the time period 
of transition between two steady states of the queue when the 
rate of the unresponsive flow changes. 

From modeling perspective, the study of the transient be- 
havior of CHOKe is an arduous task for two main reasons: 
(a) leaky nature of the queue, meaning that packets already 
in queue may be dropped later, (b) continuous state transition 
of the queue in the transient regime. Due to (a), the delay 
of a packet is not merely the backlog the packet sees upon 
arrival divided by the link capacity as in non-leaky queues. 
Besides, the spatial packet distribution of a flow in the queue 
can be nonuniform throughout the queue. Due to (b), many 
parameters that characterize the queue (e.g., flow matching 
probability, backlog size, skewed packet distributions of flows 



in the queue) are likely to be dynamically changing. Both 
constraints prohibit us from making "safe" simplifying as- 
sumptions in analyzing the transient behavior of CHOKe. 

In this paper, we take the first step in characterizing the 
transient UDP transmission rates at a CHOKe queue in the 
immediate aftermath of change in the UDP traffic arrival rate. 
In particular, we focus on: (1) how UDP utilization evolves 
in transient time, and more importantly (2) the limits, i.e., 
how far the utilization goes up or down in transient time. We 
notice that as the UDP traffic arrival rate goes up or down, 
its transient transmission rate can go in opposite direction in 
a dramatic fashion. For example, when the UDP rate sharply 
decreases, its utilization rapidly soars and often exceeds the 
steady state limits asserted above. Extreme transient behaviors 
are observed when a very high rate UDP flow abruptly stops. 
In such cases, UDP transient utilization can suddenly jump 
from 0% to over, say, 70%. This intriguing phenomenon 
cannot be explained with the literature results. 

The contribution of this paper is several-fold. First, the 
above intriguing phenomenon is illustrated with examples, 
which forms a motivation of the work. Second, for any given 
UDP traffic rate, we derive quantitatively the queue parameters 
that characterize the spatial properties of the queue in a steady 
state. Third, by leveraging the queue parameters above and 
abstracting the UDP rate change by a factor, both the evolution 
and the extreme points of UDP throughput in transient regime 
are analytically derived for any arbitrary UDP rate change. 
Last but not least, we obtain generic plots that succinctly 
represent both transient and steady state UDP throughput 
behaviors. Extensive simulations confirm the validity of the 
theoretical results. 

D. Structure of paper 

The rest of this paper is organized as follows. Next section 
presents the system setup, basic assumptions and notation. 
Sec. [3] explains the motivation using examples, and presents 
some background on steady state behavior. Since the transient 
phase is a transition period between two steady states, dis- 
cussion on steady state models is as relevant. Our theoretical 
models are presented in Sec.|4] In particular, Sec. 14- Al derives 
the insightful rate conservation argument and obtains further 
simplifying assumptions required for the transient analysis; 
Sec. 14-BI lays out the theoretical foundation on the spatial 
distribution model just before rate change; Sec. 14-CI tracks 
the UDP link utilizations and derives its properties during 
the transient period. Sec. [5] presents model validation and 
simulation results. Sec. [6] concludes the paper. 

2. System Model and Notation 

The studied system is shown in Fig. Q] where N rate- 
adaptive similar TCP flows share a link with a single un- 
responsive / aggressive UDP flow. Flows are indexed from 
0, • • • , N, where denotes the UDP flow. Since TCP flows 
are similar, hereafter 1 denotes a typical TCP flow. Packets 
are queued and scheduled in the FIFO manner on to the link. 
The steady-state backlog size in packets is denoted by b. We 
assume large N and b. 

The full notation is summarized in Table Q] The key 
performance metric is the flow (link) utilization denoted by 
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Fig. 1: System model. 



Hi, i € [0,1]. It represents the flow's share of bandwidth on 
the link. Given UDP link utilization jj.Q and link capacity C, 
the throughput of all flows can be computed. 

A queue parameter that is not indexed with a control 
variable, say time t, designates the value of the parameter in 
the steady state. Otherwise, the parameter is dependent on that 
control variable. For example, fio designates UDP utilization 
in a steady state, while /Xo(f) is changing with t, more likely 
during a transient regime. 

TABLE I: Notation 



Symbol 


Description 


N 


number of TCP flows 


r 


congestion(RED)-based dropping probability, or 
ambient drop probability (common to all flows) 


■v, 


source rate of flow i 


V-i 


link utilization of flow i 


bi 


amount of flow i packets in buffer 


b 


total backlog b = Y,f-n bi in packets 


hi 


the ratio bj/b (matching probability) 


ye [o,b] 


position in queue 


v(y) 


packet velocity at y 


<y) 


queueing delay to reach at y 


Pi{y) 


prob. of finding flow f 6 [0, 1] at position y 



A consequence of large N is that for a TCP flow, the packet 
matching in CHOKe is rare and its drop rate is mainly due to 
congestion (RED). That is, 



b 



< 



1 



:0 



bo+Nbi ~ N 
3. Motivation and Background 
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A. Motivating Examples 

To clarify our motivation, we provide example scenarios 
employing the network setup depicted in Fig. [T] There are 
N = 1 00 TCP sources and a UDP flow whose arrival rate is 
dynamically varying. The simulation parameters are described 
in full in Sec. [5] 

Example 1: The initial UDP arrival rate jto was 0.5C and 
0.25C, where C is the link capacity, for Experiment 1 and 
Experiment 2, respectively. At t = 21, the xq suddenly jumps 
by factors of 4 and 12 to 2C and 3C, respectively, and then 
returns back to 0.5C and 0.25C at t = 22. We conduct 500 
replications of the two experiments and the resulting UDP 
flow utilizations (as measured using time intervals of 10ms) 
are shown in Fig. [2] 

Example 2: The initial UDP arrival rate was IOC. As shown 
in Fig. [3] the input UDP rate flaps between 1C and 10C every 
250ms. The figure shows UDP utilization averaged over 1000 
replications, with measurements taken every 1ms. 
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Fig. 2: Transient UDP utilizations when rate xq changes by a 
factor of 4 and 12. 
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Fig. 3: Transient UDP utilizations when xq flaps between 1C 
and 10C every 250ms. 



B. Observation and Objective 

In this paper, we are interested in the peculiar behaviors 
(e.g., shaded in Fig. |2]i when UDP arrival rate varies. The 
behaviors can represent transitions of a CHOKe queue from 
one steady-state to another and typically lasts a full queueing 
delay. We note during this transient phase that, 

• The queue does not simply undergo a smooth transition 
between the two steady states. Surprisingly, the change 
in transient UDP utilization is not determined by the j^o 
values of the two steady states. Rather, it follows the 
direction opposite to the arrival rate change. Specifically, 
if UDP arrival rate goes up, its utilization dips and 
vice-versa. For example, consider Experiment 2 shown 
in Fig. |2] Instead of steadily increasing from a steady 
state value Ho.o.isc = 16% to another steady state value 
Mo.3c = 21%, the transient UDP link utilization first 
whittles down to 3.75% at t = 21.18s. (See Fig.|5]to find 
out steady state Hq values.) 

• While UDP buffer and utilization bounds stipulated in 
earlier works IIT41 . ifTBl hold in the steady-state, they are 
easily violated during the transient phase. For example, 
Mo < 26.9 for Vxo during steady-state. However, during 
the transient phase, /irj may increase abruptly by several 
factors. In Fig. [3] for instance, jUo(21.7i) =72%. 



4 



What are the lowest and highest UDP utilizations during the 
transient phase? Note that these extreme UDP utilizations dur- 
ing the transient regime depend on the measurement interval 
/ window. With large windows, the transient behavior gets 
diffused and evened out by the adjacent steady-state results. 
The objective of this work is to obtain the extreme packet- 
level UDP utilizations (i.e., the highest or lowest utilizations 
as measured for each packet) by analysis. 

These observations have crucial and practical implication. 
Internet flows, dominated by short Web transfers, often see 
fluctuating available bandwidth and may adapt their sending 
rates. For example, in periods of high UDP link utilizations, 
the number of Web document transfers may dwindle and their 
transfer completion times may get inflated. 

Since the transient regime represents a departure from a 
stable queue state, its behavior may be influenced by the earlier 
steady state, as we shall soon see. Therefore, it is fitting to 
briefly discuss the steady state behaviors. 

C. Background on Steady State Models 

As mentioned earlier, steady state CHOKe models, e.g., 
03), 03], are important for predicting transient behaviors. 
We call the model in 031 the overall loss model and the one 
in |[T4ll the spatial distribution model. Both models assume the 
CHOKe dropping and RED-based dropping are reversed for 
analytic simplicity, as shown in Fig. [4] The only independent 
parameter in both models is the UDP flow arrival rate, xq. 
Thereupon, change in xq causes the departure from current 
UDP utilization and kicks start the transient regime. 




0.1 0.25 0.5 1 2 3 5 10 
UDP arrival rate x /C 

Fig. 5: Steady-state UDP utilization jUq and buffer share Iiq. 



includes the packet velocity and the probability of finding a 
flow packet at any position in the queue. An important concept 
introduced is thinning which refers to the decaying of UDP 
packet velocity as the packet moves along towards the queue 
head. This steady state model, and the notion of thinning, are 
useful for studying the transient properties as well. We present 
the basic model, and provide pertinent assumptions that enable 
us to extend the model to the transient regime in Sec. [4] 

4. Modeling the Transient Regime 

Before delving into the transient model, we explain the 
transient behaviors using an argument, which we call the rate 
conservation law. This argument also provides an insight to 
extending the SDM to transient regimes. 
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Fig. 4: Schematic of dropping in CHOKe. 

1) Overall Loss Model (OLM): This model derives the 
steady-state UDP utilization (Iq and buffer share Iiq by first 
deriving the flow loss probability incurred by both CHOKe 
and RED parts of the CHOKe queue. The analysis results in 
a nonlinear numerical relation between ho and UDP input 
rate xq given by ©, (01, and is graphically shown in Fig. [5] 
Further details can be found in |[T5l . 



Mo 



ln[(l-fc )/(l-2/*o)] 



[(l-*o)/(l-2Ao)]+ln[(l-Ao)/(l 
xo(l-r)/C = n /(l-2h ) 



■2h )} 



(2) 



(3) 



The figure verifies the limit and asymptotic properties 
reviewed in Sec. Il-Bl The steady state simulation results in 
Fig. |2] (areas outside the shade) also match /Xo values that can 
be found from Fig. [5] 

2) Spatial Distribution Model (SDM): To characterize the 
peculiar transient properties of a CHOKe queue, Tang, Wang 
and Low developed a novel ordinary differential equation 
model fl4*l . The model captures the spatial distribution of 
flows in a CHOKe queue using queue properties at the tail 
and head as boundary conditions. This spatial distribution 



A. Rate Conservation Argument 

Before a UDP packet can be admitted into a CHOKe queue, 
it must survive both the RED and the CHOKe based dropping. 
The probability of packet admission into queue is then (1 — 
r)(l — ho(t )). Once in the queue, UDP packets can still be 
lost. This is because incoming packets that evade RED-based 
dropping (with probability 1 — r) may trigger flow matching 
(with probability ho(t)) and cause dropping of the matched 
packets. In addition, UDP packets can also leave the queue 
due to transmission with rate jUo(f)C. Summarizing, we get 
a system invariant that captures the rate of change in UDP 
buffered packets as follows: 

db (t) 



dt 



= x (t)(l-r)(l-h (t))-x (t)(l-r)h (t)-^(t)C 
= Xb(t)(l-r)(l-2Ao(0)-/io(t)C (4) 

Let us call © the rate conservation law. That is, the 
rate of change in UDP buffer occupancy is the difference 
between the flow queueing rate xo(?)(l — — ho(t)) and 
the outgoing rate. The outgoing rate in turn is the sum 
of the departure/transmission rate given by jio(t)C and the 
leaking rate given by xo(f)(l —r)ho(t). Here the leaking rate 
denotes the rate with which a queued UDP packet matches the 
incoming packet and is consequently dropped. 

The corresponding equation for a TCP flow is, where note 
hi = due to large number assumption (see (Q])), 

db\{t) 



dt 



= xi(t)(l-r)-Hi(t)C, 



(5) 
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where, since the link is fully used, 

l-Mo(0 



Mi(0 

Since, trivially, b(t) = b (t) • 



N 

-Nbi(t), 



we get 



db(t) db (t) „db\(t) 



N- 



(6) 



(7) 



dt dt dt 

We remark that during a stable / steady state db((t)/dt « 
0, ie (0,1) (see Fig.©. 

Now, let us assume an abrupt change in UDP arrival rate xq 
and note the following in the immediate aftermath: TCP flows 
react slowly in response to the sudden change in UDP arrival 
rate. Notably, TCP flows react in an RTT timescale, but the 
transient behavior lasts for sub-RTT scales. Hence, during the 
transient phase, 

db\(t) n db{t) db (t) 



dt 



:0 



dt 



dt 



(8) 
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(a) Experiment 1 (0.5C/2C). 
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(b) Experiment 2 (0.25C/3C). 

Fig. 6: Transient regime: db/dt « dbo/dt and db\/dt « 0. 
Steady state: db/dt « 0, dbo/dt « 0, db\/dt « 0. 

Fig-Elplots the simulation results for the two experiments of 
Sec. l3-AI The figures confirm the argument above and Eq. ([8]). 
As can be seen, for most parts of the simulations, the rate of 
change in total buffer occupancy b is almost solely due to 
change in UDP buffer occupancy bo. 

We warn the reader that <[8j may not hold outside the 
transient regime. For instance, when the UDP flow arrival rate 
plummets at t = 22s, TCP flows respond by increasing their 
sending rates after around a round-trip delay (see especially 
Fig. |6b]i. Subsequently, dbi/dt ^ 0, but dbo/dt « 0. From ©, 
db/dt K.Ndb\/dt. That is, db/dt swings from dbo/dt during 
the transient phase to Ndb\/dt following the transient phase. 
However, since the TCP flows are largely in a congestion 
avoidance phase, the rise in Ndb\/dt is not as significant as 
that of dbo/dt in the transient phase. After absorbing TCP 
bursts for a while (in a few round-trip cycles), the queue 
eventually settles to a new steady state determined by the new 
UDP arrival rate, and then dbi/dt w0, ; G (0, 1) once again. 



Remark: Combining © and ©, we conclude that TCP 
packet arrival rate to the CHOKe queue matches its transmis- 
sion rate during the transient regime. 

Now we are in a position to explain the transient behaviors 
shaded in Fig. [2] and Fig. [3] from the perspective of the rate 
conservation law. Rearranging (0), we get 

*>(')(! -')„ 1 db o(' 

C 



■(l-2Ao(0)- 



(9) 



C v y " C dt 

Interestingly, (O captures all pertinent behaviors of the 
system. In the steady state, dbo(t)/dt « 0, and 



Mo = 



*o(l-r) 
C 



(1-2/zo) 



(10) 



which is ([3). dTOb is a key equation in the OLM model and 
can reproduce the steady state results depicted in Fig. 

For the transient behaviors, we explain only the dips in jjq 
shown in Fig. [2] but similar arguments follow for the peaks as 
well. An abrupt injection of UDP rate xq at t = 21s rapidly 
ramps up the second term on the right hand side (r.h.s.) of 
(O (see also Fig. [§}. Despite the rapid rise of xq(1 — r)/C, its 
contribution to jllo(f ) is counteracted by a corresponding rise of 
h (t) = b Q {t)/b(t) (see Fig. |7] where h (t) ->45%). Therefore, 
IM)(t) is mainly influenced inversely by the rate dbo(t)/dt (cf. 
Fig. [6] and Fig. [2j. Specifically, focusing on Fig. [2] and using 
measurement intervals of 10ms, we observed the following 
: when xo increases 12-fold at t = 21s, dbo/dt — s- 3500 (as 
shown in Fig. [6bb but the utilization jUo — > 3.75%. Conversely, 
when xo slackens by a factor of 12 at t = 22s, the dbo/dt 
falls to 1000 below (as shown in Fig. [6bJ but the UDP link 
bandwidth share soars to 56.5%. 




20.5 



21.5 
Time (s) 
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Fig. 7: At t = 21, UDP buffer share ho jumps radically, 
nullifying the impact on jUo of the sudden change in xq. 

While the rate conservation argument explains the transient 
queue dynamics well, due to several dynamically changing 
parameters (/zo(f), bo(t), b{t)) in the transient regime, it is 
difficult to derive quantitative jj$ results directly from (|9). 
Nevertheless, some of the insights we gained become useful 
for extending the SDM model to the transient regime, as we 
shall soon see. 

B. Modified Spatial Distribution Model 

This section is dedicated to the development of the SDM 
in concert with the transient regime. A schematic diagram of 
this model is illustrated in Fig. [8] 
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Fig. 8: Schematic diagram: Decay of UDP velocity Po{y)v(y) 
in queue. Look at the similarity to Fig. [9b] with xq = IOC. 



1 ) Model parameters: SDM can be described by a few key 
parameters (see Fig. [8] and Table H). The parameters are the 
queue position/point/slot y 6 [0,b], the packet velocity v(y) at 
y, the probability p,(y) of finding a flow i packet at y, and the 
queueing delay x(y) for the packet at the tail to arrive at slot 
y. Queue position y is indexed from tail to head as {0, • • ■ , b}. 
The packet velocity v(y) is the speed with which packets move 
towards the head of the queue, and is defined as: 



v(y) = dy/dt. 



(11) 



The packet velocity at queue tail v(0) is simply the full 
queueing rate J^LqX^I — r)(\ — hi) (see Fig. |4). At queue 
head, however, v(y) is merely the link capacity, i.e., v(b) = C. 
The packet velocity v(y) is related to the queueing delay x(y) 
accumulated in going from tail y = to slot y as follows, 



dt = dy/v(y) 



J v(s) 



(12) 



where the equation on the left side is obtained from (fTTT i. 

Of course, the full queueing delay is %{b). Alternatively, it 
can be derived using queueing principles. The rate of departure 
of TCP packets is C(l — jUo), and average number of TCP 
packets in queue is given by b(l — ho)- As far as the TCP 
flows are concerned, the model is a non-leaky queue (since 
h\ = 0). Therefore, we apply Little's law to obtain, 



x(b) 



b(l-hp) 
C(l-Mo)" 



(13) 



Another useful spatial parameter is p;(y) — the probability 
of finding a flow ; packet at slot y. Trivially, 



Po(y)+N Pl (y) = l ye [0,1 



(14) 



Here, piiy) is closely related to the packet velocity v(y): It 
quantifies the fraction of flow i's packet velocity at y to the 
total packet velocity v(y). For instance, at queue head y = b, 
Pi(b) = fii, i.e., the probability is simply the flow utilization. 

Summarizing the two important parameters pi(y) and v(y) at 
queue tail and head, the following boundary conditions apply, 
which are also illustrated in Fig. [8] Here, we ignore TCP flow 
matching as discussed earlier for ([TJi, i.e., h\ w 0. 



v(0) =*<,(! -r)(l-*o) +Nx l (l-r)(l-h l ) 

wjco(l-r)(l-*o)+iVxi(l-r) (15) 
x {l - h )(l - r) x Q (l-ho) 



Po(0) = 



(16) 



P o(b)C 
= PoC 

N Pl (b)C 
= (l-Ato)C 



[x (1 - ho) +Nxi] (1 - r) xo(l — ho) + Nxt 
Pi(0) = — *\ , 



x (l-h ) 
v(b)=C, p (b)=Ho, 



Pi(fc) = Mi = 



1 -juq 
N 



(18) 



Note that pi (b) in ( fT8l is the same as ([6). As noted in 
Sec. 14- Al TCP transmission rates do not change during the 
transient regime, resulting in constant TCP packet velocities 
throughout the queue. In a congested CHOKe queue with 
high UDP rate xq, however, the total packet velocity v(y) is 
continuously decreasing because UDP arrivals trigger packet 
drops through flow matching. In fluid terms, we say the UDP 
fluid gets thinned as it moves along the queue. We formalize 
the notion of thinning and use it to derive the slot parameters 
next. 

2) Ordinary differential equation model: The UDP portion 
of the packet velocity at the tail (y = 0) is given by po(0)v(0) = 
xo(l — /io)(l — r) and the amount of UDP fluid in small time 
dt at the tail by po(0)v(0)iif . The corresponding values at y 
are Po{y)v(y) and po(y)v(y)dt, respectively. Traveling from 
queue tail to y takes t(y) during which time xo(l —r)x(y) 
new packets would arrive to the queue. Each arrival triggers a 
flow matching trial and drops the small volume of fluid with 
success probability 1 jb. The probability that the UDP volume 
escapes matchings by all arrivals is (1 — l/t>) x o( l ^ r ) r (y) . The 
UDP packet velocity at y is therefore thinned or weakened as, 

Po(y)v(y) = Po(0)v(0)(l - i/6)*(i-')*M (19) 

On the other hand, there is no thinning for TCP and hence a 
constant TCP packet velocity throughout. 

Npi (y)v(y) = Npi (0)v(0) = (1 - p )C (20) 

During the transient regime, ( f20b is still valid due to the slow 
reaction of TCP congestion control, as discussed in Sec. 14- Al 
Rearranging ( f20b and using dT4b . we get 



v(y) 



Pi(0)v(0) (l-Mo)C (l-/io)C 



(21) 



Pi 00 Npi(y) 1-poW 

Define parameters a and j3 as follows, which will often be 
used throughout the rest of the paper: 

_ l-po(O) 
°' Po(0) 
J3 : = ln(l-l/Z>) 

Using ( fT4b . ( fT3] l, ( fT6] l and ( |20l i, a can equivalently be 
rewritten as, 

1 - Po(0) = Npi (0) = Npi (0) v(0) = (l-po)C 
p (0) p (0) p (0) v(0) xb(l-r)(l-Ao) - 
Taking logarithm of key equation (T% first and then differ- 
entiation w.r.t. y, we get 

ln(p (y)v(y)) = ln(p (0)v(0)) +x (l - r)j3r(y) 



Po(y) , v Jy) 

Po(y) v(y) 



= 0- 



xo(l-r)j3 

v(y) 



(22) 



7 



where t' (y) = l /v(y) from (fT2l i. 

Applying ( f2TT > for v(y) and v (y) and inserting into l.h.s. 
of d22l) . the following ordinary differential equation (ODE) is 
obtained. 



PoM . Pol^) _xo(l— r)j8 



(23) 



PoM 1-P0M y M 

Eq. d23l establishes a foundation for the analysis in the 
remaining part of the paper. 

First, solving PoM from ODE (l23l (see Appendix for the 
proof), we have the following relation between PoM and i(y): 

Lemma 1. 



PoM = 



e x Q (l-r)PT(y) 
, +e x (l-r)pt(y) ■ 



(24) 



In addition, substituting (|2TT > for v(y) in (l23l gives, 
PoM 1 PoM *o(l-r)fl / \\ 



Dividing by (1 — PoM)> we obtain 
PoM PoM 



_ xo(l-r)/3 

PoM(i-poM) ' (1-P0M) 2 (i-Mo)c 



Upon integrating w.r.t. y, we get 



(l-Mo)C 

Jro(l-r)j3 
(l-Mo)C' 



PoCO 





In 



Po(s)(l-Po(s)) 
PoO) 



PoM 



(1-PoW) 



1 " PoO) 



Let K = ^-M. We then get: 



Ky = ln 



PoM i-po(o) 



1-poM Po(0) 



1 - PoM 



1 



1 



1-poM i-Po(o) 



(25) 



from which, the following relation between y and p {y) is 
easily established: 



Lemma 2. 

1 



K 



■In 



qpoM 
1 - Po(.v) 



PoM-Po(O) 



tf(l-poM)(l-Po(0)) 



(26) 



where K = and as above, a := - ( l -^ c 



(l-Mo)C " J "^ v " • Po(0) ^o(l-r)(l-ft ) ■ 

Remark: Lemma [TJ and Lemma [2] are crucial, through 
which, the slot and its associated parameters are inter-related 
as defined by p (y) ([24), y ©, vM <ED and r(y) OH). 
However, from the model in lfT4l . it is difficult to obtain for 
each slot y corresponding explicit expressions of the spatial 
parameters PoM> v(y) and T(y), hence limiting the application 
of results in |[T4l to the transient analysis. This is because 
PoM an d v (y) were given in terms of queueing delay t(y), 
but no explicit relations were afforded between queue slot y 
and anyone of the triplets (v'M,poM, I'M) m El- 
Remark: More importantly, the analytical principles lead- 
ing to Lemma [TJ and Lemma |2] also apply to the transient 
regime, making the later transient analysis possible. 



3) Properties of queue dynamics: This section briefly out- 
lines properties of queue dynamics PoM, v(y) and l(y). The 
detailed proofs are in the Appendix. 

We start with PoM- To prove the properties of PoM in 
Lemma 21 we need the following intermediate result. 



Lemma 3. 



l<i^2<2 
" I-/10 " 



The following properties hold for PoM> an d f°ll° w from 
( f24b and Lemma [3] 

Lemma 4. Given xq > 0, 

(a) If xo ^ C/2, PoM i s strictly convex decreasing. 

(b) Otherwise, PoM is concave decreasing from po(0) to p ( j, 
and convex decreasing from p^ to po(^) where the critical 
point Mi Po) ^ s gi yen by 



Po* =PoM) = 3 



II - 

k n \2j 1 a: 2(1 



1 l-3po(0) 



Po(0))' 



For I'M' a PPlyi n g d24b to (fJJJi first, then the following 
property can be verified from (|2"TT i 

Lemma 5. Given xq > 0, and boundary values v(0) and v{b) 
defined respectively by $15[ and (fT8l l, the packet velocity v^) 
is convex decreasing throughout the queue. 

For x(y), the following property holds. 

Lemma 6. Given xo > 0, t(0) = and z{b) defined by ( fT3l , 
queueing delay %{y) is strictly convex increasing. 

To visualize the queue dynamics, Fig.|9]plots for each queue 
position y the probabilities po(y), the packet velocities v(y) and 
the queueing delays t(y). For low and moderate rate xq, the 
spatial properties v(y) and poM 1°°^ uniform in queue, and 
the queueing delay t(y) is approximately linearly rising, like 
in a regular non-leaky queue. With increasing rate xq, however, 
the spatial distribution becomes increasingly asymmetrical in 
queue, with most of the UDP packets piled up closer to the tail 
and the packet velocity sharply decreasing to C as we move 
towards the head of queue. In such cases, since the packet 
velocities are nonuniform in queue, the queueing delays are 
also nonuniform. When xo — > °°, the queueing delay becomes 
5^ — half of the queueing delay possible in a non-leaky queue 
with the same backlog size b. 



C. Analysis on the Transient Behavior 

The transient behavior is a transition between two stable 
queue states excited by a change in UDP arrival rate. Without 
loss of generality, let us assume that the UDP arrival rate 
changes from xq to xq2 at time t = 0. At that instant, the 
snapshot of the queue exhibits the steady state characteristics 
defined by xo (see Fig. [8). We are interested in transient 
behaviors that surface within a time of %(b) (the full queueing 
delay) after rate change. Thereafter, we assume the queue 
enters the steady state defined by new rate xo2- We shall focus 
on the transition regime in [0, t(&)]. 
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x„ = 0.25C 




x a = 2C 




x a = IOC 







(a) UDP packet probability vs queue slot. 



i 0.26 0.46 0.64 0.86 
Queue slot y 

(b) Packet velocity at queue slot. 




(c) Queuing delay to reach at a queue slot. 



Fig. 9: Spatial characteristics of a CHOKe queue under different intensities of input UDP arrival rates. 



While Lemma [TJ and Lemma [2] lay a foundation for the 
analysis, we still need to make assumptions to realize it. 
Specifically, we boldly assume, during [0, x(b)], 

I: the buffer occupancy remains constant at b, and 
II: TCP arrival rates remain the same. 

Remark: It is worth highlighting that while Assumption ITT1 
is valid or accurate based on our observation of slow TCP 
reactions in Sec. 14- Al Assumption U may be very conservative 
(see Fig. [7]). Incidentally, this assumption may be the cause 
of approximation errors in the analysis. Nevertheless, results 
based on these assumptions are highly satisfactory. 

We are now ready to conduct the analysis. Consider at time 
t = the UDP fluid at position y. It has a velocity po(y)v(y), 
and a remaining age in queue of x(b) — x(y). Owing to FIFO, 
the UDP flow transmission rate at time x(b) — x(y) is due 
to the transmission of this UDP fluid. By transmission time 
x(b) — x(y), the UDP velocity has been thinned according to 
the new rate X02 as: 



With Assumption [II] the total packet velocity v at time 
x(b) — x(y) is the sum: 

v(x(b) - x(y)) = p (y)v(y)(l - l/^M^M) 

+ (l-Mo)C (29) 

where the second term represents the velocity contributed by 
the TCP flows, given by d20l i. 

The instantaneous UDP link utilization po follows from i 
and ( |29l simply as, 



po T (b) - T (y = T(h) I ^ (30) 

v(x(b) - x{y)) 

Note that, with 1211 . which holds for both the steady-state 
and the transient regime, we can express the UDP packet 
velocity, p (y)v(y), as 

Po(y) 



pQ{yHy) -l- P0 (y) 
with which, we further obtain 



(l-Mo)C 



(31) 



Po(y)v(y)(l - i/ b y^-^)-r(y))_ (27) ^ m T(y)) = 



1-PoW 



1-1 

b 



-■*m(T(A)-T(y))- 



Remark: d27l i is similar to (fT9l . Unlike (fT9l , however, there 
are two distinct parts of thinning for the UDP fluid located 
at slot y € [Q,b] at t = 0. First, in going from tail to y, the 
thinning is according to Xq, and the duration of thinning is the 
queueing delay so far, i.e., x(y). This part is reflected in (|27] | 
by Po(y)v(y). Second, for the remaining duration x(b) — x(y), 
the thinning is due to the new rate xo2- Overall thinning, before 
eventual transmission, of the UDP fluid found at slot y at t = 
is then, 

po(0)v(0)(l -l/b) Xo{l - r)T{y) (l - i/j,)*te( 1 -' , )(*( 6 )-*G0). 

Above, we use the same RED / congestion-based drop 
probability r even when the UDP arrival rate changes. In fact, 
our extensive simulations show that r is often insignificant and 
can be ignored altogether. CHOKe's excessive flow dropping 
keeps the average queue size avg in check, and this in turn 
lowers the r in comparison to that in plain RED. Hereafter we 
choose to ignore r. That is, 



po(y) 



-i 



1 + ae -- r 02/3T(^)+/3T(v)(j:o2-^o) 



-1 



(32a) 



(32b) 



In obtaining the reduced forms (|32at , (|32bt , we used (|3TT i and 
(l24l respectively. 

d32bt captures the evolution of UDP utilization during the 
transient regime. We summarize it in the following lemma. 

Lemma 7. Assume at t = 0, UDP arrival rate changes from 
to xq2- The UDP link utilization at time AT G [0,x(b)} is 
given by, 

-1 



Ho(AT) = 



1 



-x ,Jz{b)+P(z(b)-KT){x 02 -x a ) 



r «0. 



(28) 



Here, we would like to remark that similar observation or 
assumption has been made for CHOKe analysis in the litera- 
ture IfBI . 



where, b is backlog size at t = and x(b) is given by ( fT3l . 

It is trivial to see that po(0) = Po- In addition, it is easy to 
prove that when x 02 = x , po(Ar) = [i for VAT G [0,x(b)]. 
The proof is similar to the proof of Lemma [8] 

Theorem 1. Assume a CHOKe queue characterized by steady 
state UDP probabilities po(y), y G [0,b] and input UDP rate 
xo- Further assume the UDP rate changes to xq 2 > at t = 0. 
(a) The transient UDP utilization is upper bounded by Po(0). 
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(b) This upper bound can be achieved when xq2 — and at 
time x(b). 

Proof: Since x(b) ^ r(y) (see Lemma [6]i and /3 < in 
(I32al i. we have exp(— XQ2^('c(b) — t(v))) ^ exp(O) = 1 and 



fi (T{b) - z{y)) < 



1 



1-PoW 



PoGO 
<Po(y)<po(0) 



(33) 



In (l33l l. we used the property that po{y) is a decreasing 
function (see Lemma |4|i to state that /Xo(t(£>) — T(y)) ^ Po(0). 
When x Q2 = in (|32a>, ;Uo(tO) - r(y)) = p (y). See Fig. [TO] 
for the relationship between Po(y) in queue and transient UDP 
rate. ■ 



M0X&))=P0(0) 




t = 



Fig. 10: Relationship between steady-state po(y) and transient 
UDP utilization /io(f ) when the flow stops, i.e., xq2 = 0. 

Theorem Q] states that if the UDP flow stops (xq2 = 0), it 
will attain exactly the utilizations po(y) in reverse order of 
time shown in Fig. [8] That means, the steady-state probabilities 
Po(y) associated with the previous UDP input rate xq succes- 
sively turn out as transient utilizations when the flow stops. 
Therefore, the transient utilizations increase from po(b) = jj$ 
at t = to p (0) at t = x{b). 

Theorem Q] vindicates our choice to leverage steady-state 
CHOKe model for understanding the transient behavior. 

So far we have discussed the two special cases when xq2 = 
xq (no change) and X02 — 0. Now consider a scenario where 
xo2 ^ {0,xo}, in particular X02 —> °°. It is plausible that when 
XQ2 is high, the exponential terms in d32ab . or equivalently 
in Lemma [7] may become so large that the utilization may 
quickly plunge to very low values. From our observation, we 
remark that the UDP utilization in the transient phase moves 
in the opposite direction to the change of UDP input rate that 
triggers the phase. Since the UDP rate change impacts the 
rate dbo/dt, the above observation is coincidentally the same 
as the one noted in Sec. 14- Al 

We now generalize our findings and obtain the extreme 
values. Note that from Lemma |7J [Iq(AT) is decreasing or 
increasing with AT, depending on whether X02 is greater than 
xq or not, so the extreme is obtained when AT = i(b). In other 
words, the extreme (lowest or largest) values occur i(b) after 
rate change. That means, the last packet of the old rate xq is 
transmitted with the extreme utilization jj.^ = /Io(t(£>)). The 
next theorem gives the maximum or minimum value. First, 
we start with Lemma [8] which captures the special case when 

05 = 1. 



Lemma 8. For a = 1, Ho(i:(b)) = Ho- 

Theorem 2. Assume the current steady-state utilization /Xo 
when the input UDP rate is xq. If X02 = 05*o at / = 0, Ct€ [0, °°), 
then the extreme (minimum / maximum) UDP utilization 
during the transient regime is given by 



1 +a 



1 - Mo 
alio 



-l 



(34) 



See the Appendix for proofs of Theorem [2] and Lemma [8] 
Corollary 1. For a = 0, or when the UDP flow stops, 

M0(T(fo))=p (0). 

Theorem [2] is a key finding of this paper and can be visual- 
ized using Fig. QT| Note that the figure illustrates both steady- 
state and (extreme) transient UDP utilizations for the selected 
UDP arrival rates. When a = 1, by Lemma [8] the values 
shown are the steady-state utilization /Iq for the given UDP 
arrival rate xq. When the UDP arrival rate xq changes (i.e., 
a ^ 1), the graph shows how far the UDP utilization can go 
up/down when xq abruptly decreases/increases, respectively, 
to xq2- For example, assume an initial UDP rate of xq = 2C. 
The utilization for this input is read from the figure at a = 1 
(also from Fig. 13) as /i = 25.0%. When x Q2 = 0.2C which 
corresponds to rate change by a factor of a = 0.1, the transient 
utilization surges to jio(T(b)) =56.5%. If the flow stops, by 
Corollary [T] the utilization instead jumps to Po(0) =60% (see 
also Fig. [9aY Similarly, when a flow of initial UDP arrival 
rate xq = 3C stops, the transient utilizations can surge to a 
whopping 67% from the initial 21%. 



= 3 



a. 
Q 




a 



X02/X0 



Fig. 11: The impact on extreme UDP utilization of rate factor 
a. Five previous inputs x E {0.25C,0.5C, 1C,2C,3C} are 
shown. 



5. Evaluation 

In this section, we validate the results using simulations 
performed in ns-2.34. The network setup shown in Fig. Q] with 
the following settings is used: C = 20Mbps or 2500 pkt/sec, 
link latency Ims, buffer size 1000 packets, N = 100 TCP 
flows each of type SACK, RED buffer thresholds (in packets) 
min r /, = 20 and max^ = 1000. Packet sizes are 1000 bytes. 
Flows start randomly on the interval [0,2] sec. 
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We conducted extensive experiments, each simulation repli- 
cated 500 times if not otherwise highlighted. The 95% con- 
fidence intervals are so small that they are not reported. 
We remark that in computing the simulation results, unless 
otherwise stated, we have used a time window of 1ms. Since 
C=2500 pkts/sec, this is 2.5 times more than the per packet 
transmission interval assumed by the model, but the error due 
to this disparity is small and can be ignored. Section 15-AI 
presents the validation of the model, and Section |5^Bl presents 
additional simulation results. 



A. Model Validation 

In this section, we validate the two important results of this 
paper: Theorem [2] and Lemma UJ 

1) Validation of Theorem |2J In Fig. [T2J we show for 
selected initial UDP arrival rates the impact of rate change 
by factor a e [0.01,10]. As can be seen, the simulation 
results accurately match the model predictions. For instance, 
for xq = 3C and X02 = 0.03C, or a = 0.01, the maximum 
utilizations obtained by the model and simulations are 67% 
and 65%, respectively. 



Q 

3 



UJ 




Fi 



a = x 02 /x 

g. 12: Validation of extreme utilization stated by Th. [2] 



2) Validation of Lemma [7} For rate factors of a = 
5, 1/5, 10, 1/10, Fig. [131 shows the evolution of transient 
UDP utilizations obtained through simulation and the analyt- 
ical model stipulated by Lemma |7j The steady state backlog 
size b required for the theoretical plot is taken from the steady 
state simulation just before rate change. 

As discussed earlier, there may be two sources of approxi- 
mation errors for the theoretical results. First, the conservative 
assumption of constant backlog size during the transient phase. 
Second, the difference in measurement intervals. The model 
results are tallied per packet transmission time (0.4ms), while 
simulation results are based on 1ms interval. Despite these 
differences, the model and simulation results are reasonably 
matching, even more so for moderate initial xo (see Fig. 1 1 3ab . 
In all figures, the approximation errors are negligible at the 
beginning. As we move further in time during the transient 
phase, however, the assumption of constant buffer occupancy 
b fails to hold. Consequently, the errors become larger. In the 
figure, the model can be treated as a bound of the respective 



0.25 

=P 0.2 
o 

■§ 0.15 

_N 

1 0.1 
Q_ 

§ 0.05 



20.95 



Simulation 

— Model 
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Time (s) 

(a) xq = 0.25C, a = 5 



21.15 



21.2 




20.95 21 21.05 21.1 21.15 21.2 

Time (s) 

(b) xq = 3C, a = 1 /5 



Simulation 

— Model 




20.95 21 21.05 21.1 

Time (s) 

(c) x = 0.25C, a = 10 



21.15 



21.2 




21.05 21.1 
Time (s) 



(d)x = 3C, a = 1/10 
Fig. 13: Validating the transient utilization Eq. (I32bb 



transient utilization. Additionally, the extreme utilizations dur- 
ing the transient regime (i.e., the lowest utilization in Figs. I Dal 
and !13cl and highest utilization in Figs. I13bl and !13d1 ). which 
can also be verified from the generic utilization curve shown 
in Fig. [12] are very close to the simulation results. 

B. Miscellaneous Results 

The next three subsections present additional results. The 
first two return to the motivational examples discussed in 
Sec. [3] The results in Sec. l5-B3"1 are based on a different traffic 
model — Web Traffic. 

1) Results on Example 1: For the two experiments of 
Example 1 discussed in Sec. 13-AI Table [n] tabulates three sets 
of extreme UDP utilization values: theoretical values based 
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on Theorem |2j and two sets of simulation values averaged 
over measurement windows of 0.4ms and 10ms. Note that 
since transient utilization is continuously changing, the time 
granularity of measurement windows are critical (see also 
Sec. I3-Bb . In Fig. [2] the curves are based on a 10ms window. 
Compared with simulation results, the theoretical values in 
Table HI1 seem to represent lower utilization bound when a > 1 
and upper utilization bounds when a < 1. Nevertheless, the 
theoretical and simulation values based on packet transmission 
time window are remarkably close, and the bounds are tight. 

TABLE II: Lowest and highest transient UDP utilizations for 
the experiments of Example 1 in Sec. [3] 



Scenarios 


Extreme UDP utilization 


x Q /C X02/C a 


AiI A , Simulation 
Model w=()4ms w=10ms 


0.50 2.0 4 
2.0 0.50 1/4 
0.25 3.0 12 
3.0 0.25 1/12 


6.7% 7.5% 10.5% 
50.8% 49.3% 44.5% 

1.3% 2.0%) 3.8% 
63.2% 61.4% 56.7%) 



2) Results on Example 2: Next, for the motivating example 
in Sec. 13-Al where the UDP arrival rate alternates between 1C 
and 10C, the two extreme utilization values using Theorem [2] 
are found to be 0.015% (for change from 1C to 10C) and 
75% (for change from 10C to 1C). Similar to what has been 
observed in Example 1, the extreme values are tight bounds. 
Even using a gross measurement window of 10ms for Fig. [3] 
we still observe a peak utilization of 72%. Besides the extreme 
values, the model allows us to explore the average utilization 
in a time period within the transient regime. To demonstrate 
this, Table [HI] shows average utilizations (per every 500ms) 
using four methods: (1) steady state utilization corresponding 
to average UDP arrival rate 5.5C, (2) average of steady state 
utilizations corresponding to UDP arrival rates 1C and 10C, (3) 
average utilization based on the model, i.e., Lemma|7] and (4) 
simulation results, using measurement time windows of 1ms. 
Steady state values for the first two methods are based on the 
OLM model described in Sec. I3-C1I For using Lemma |7J the 
values for backlog size b at each 250ms interval are required. 
Our extensive simulations show that the backlog size is wildly 
changing over time. We took the backlog size at the onset of 
rate change and fed it as the input bo into Lemma [7] For 
example, at t = 21 bo = 765 and at t = 21 .25 bo = 675. As the 
table shows, the steady-state analytical bounds are far from 
the simulation results. On the contrary, the transient analysis 
nicely represents the picture of the queue even under radically 
changing traffic conditions. 

TABLE III: Average UDP utilizations for Example 2 in Sec. [3] 



Time 




Average UDP utilization 


in % 


Interval 


M0.5.5C 


Mo.ic+Mo.ioc 


Model 


Simulation 


[21,21.5) 


11.7 


14.8 


17.6 


18.9 


[21.5,22) 


11.7 


14.8 


19.1 


19.4 


[22,22.5) 


11.7 


14.8 


17.6 


18.9 


[22.5,23) 


11.7 


14.8 


18.3 


19.4 


[23,23.5) 


11.7 


14.8 


18.0 


19.3 



3) Results using Web traffic: Since the Internet flow dynam- 
ics is heavily shaped by short Web transfers, we conducted a 
500-replicated experiment explained below. The UDP arrival 



pattern is the same as in Fig.[3]except that xq = 10C for t < 21 
and xq = 1C for t > 23. The Web traffic is modeled as follows: 
Starting from t = 20s, each of the 100 TCP sources (see Fig.Q]) 
generates a Poisson process with an average arrival rate of 25 
Hz. The size of each session (file) is Pareto-distributed with 
average size of lOkB (about 10 packets) and a shape parameter 
of 1:3. This model captures the heavy tailed nature of Web file 
sizes and their transmission times [171 . Simulation lasts for 25 
seconds and over 9000 Web sessions have been generated. The 
result is illustrated in Fig. [14] 
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Fig. 14: UDP utilization in the presence of Web flows. 

Due to the huge and highly bursty Web traffic generated, 
the buffer is always full. Like in the long-lived TCP scenarios, 
UDP exhibits widely fluctuating throughput patterns during 
transient regimes. However, the extreme points of transient 
regime are generally lower in value. For example, UDP 
utilization can get as low as 0.08%, and as high as 56%. 
Nevertheless, in both cases, they are bounded, albeit loosely by 
the analytical extreme values 0.015% and 72% as discussed in 
Sec. I5-B2I We believe that the smaller values are due to higher 
ambient drop rates r (cf. r in ( f28l ) caused by RED dealing 
with persistent full buffer occupancy. To close the difference, 
further study may be conducted where the rate conservation 
argument could be exploited. 

The results in this section show that the analytical results 
and observations made in this paper may apply to a wider 
context than studied here. 

6. Conclusion 

While existing works on CHOKe reveal interesting struc- 
tural, asymptotic and limit behaviors of the queue, their 
results are limited to the steady state when the queue reaches 
equilibrium in the presence of many long-lived TCP flows and 
constant rate UDP flows. Unfortunately, they lack showing 
properties of the queue in a possibly more realistic network 
setting where the exogenous rates of unresponsive flows may 
be dynamically changing, and consequently the model param- 
eters, rather than being static, may be continuously evolving. 

This paper provides the first study on CHOKe behavior 
in the aftermath of rate changes in UDP traffic arrival. In 
particular, we are concerned with CHOKe queue behaviors 
during the transient regime which we model as a transi- 
tion from one steady queue state to another. We found that 
the performance limits stipulated in the steady state rarely 
hold for such transient regimes. Depending on the nature of 
rate change, the queue exhibits instant fluctuations of UDP 
bandwidth sharing in reverse direction. This behavior has 
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ramifications on the smooth operation of the Internet where 
most flows are rate-adaptive. Such flows may see fluctuating 
available link bandwidth and degrade in performance. By 
extending and leveraging the spatial distribution model, this 
paper analytically (1) determines the extreme points of UDP 
utilization (observed within an order of queueing delay after 
rate change), and (2) tracks the evolution of the transient 
UDP utilization following rate change. In addition, the model 
allows us to obtain generic UDP utilization plots that help 
explain both the transient extreme characteristics and steady- 
state characteristics. The analytic results have been rigorously 
validated through extensive simulations. We believe, the ana- 
lytical approach used in this paper also sheds light on studying 
transient behaviors of other leaky queues. 



Appendix 

This appendix is devoted to proving the lemmas and theo- 
rems in Sec. 14- Bl 

We need the following intermediate results. 
From (fT2l and (BTT i. respectively, we get 

T(y) = l/v(y) (35) 
v'(y)/v(y)=p' Q (y)/(l-p (y)). (36) 

Using Po(y) from (l36l l into (l22l and solving for v (y), 

v'(y)=xo(l-r)Ppo(y). (37) 

Proof of Lemma Q] 



References 

[1] J. H. Saltzer, D. P. Reed, and D. D. Clark, "End-to-end arguments in 
system design." ACM Transactions on Computer Systems, vol. 2, 1984. 

[2] M. Altaian, V. Paxson, and E. Blanton, "TCP Congestion Control," RFC 
5681, Sep 2009. 

[3] A. Demers, S. Keshav, and S. Shenker, "Analysis and simulation of a 
fair queueing algorithm," in Proc. of ACM/SIGCOMM, 1989. 

[4] S. Floyd and V. Jacobson, "Random early detection gateways for 
congestion avoidance," IEEE/ACM ToN, vol. 1, no. 4, 1993. 

[5] S. Floyd, R. Gummadi, and S. Shenker, "Adaptive RED: An algorithm 
for increasing the robustness of RED's Active Queue Management," 
Tech. Report, 2001. 

[6] S. Keshav, An engineering approach to computer networking: ATM 
networks, the Internet, and the telephone network. Addison- Wesley 
Longman Publishing Co, 1997. 

[7] A. Kortebi, L. Muscariello, S. Oueslati, and J. Roberts, "Evaluating the 
number of active flows in a scheduler realizing fair statistical bandwidth 
sharing," in Proc. of ACM/SIGMETRICS, 2005. 

[8] B. Suter, T. Lakshman, D. Stiliadis, and A. Choudhury, "Buffer Man- 
agement Schemes for Supporting TCP in Gigabit Routers with Per-flow 
Queueing," IEEE Journal on Selected Areas in Communications, vol. 17, 
no. 6, Jun 1999. 

[9] A. Eshete and Y. Jiang, "Approximate Fairness Through Limited Flow 

List," International Teletraffic Congress, Sep 2011. 
[10] C. Hu, Y. Tang, X. Chen, and B. Liu;, "Per-Flow Queueing by Dynamic 

Queue Sharing," IEEE INFOCOM, 2007. 
[11] D. Lin and R. Morris, "Dynamics of random early detection," in ACM 

SIGCOMM Computer Communication Review, vol. 27, no. 4, 1997. 
[12] R. Mahajan, S. Floyd, and D. Wetherall, "Controlling high-bandwidth 

flows at the congested router," in IEEE ICNP, 2001. 
[13] R. Pan, B. Prabhakar, and K. Psounis, "CHOKe - a stateless active queue 

management scheme for approximating fair bandwidth allocation." in 

IEEE Infocom, 2000. 
[14] A. Tang, J. Wang, and S. H. Low, "Understanding CHOKe: Throughput 

and Spatial Characteristics," IEEE/ACM ToN, vol. 12, no. 4, 2004. 
[15] J. Wang, A. Tang, and S. H. Low, "Maximum and asymptotic udp 

throughput under CHOKe," in ACM SIGMETRICS, 2003. 
[16] R. Pan, C. Nair, B. Yang, and B. Prabhakar, "Packet dropping schemes, 

some examples and analysis." in Allerton Conference on Communica- 
tion, Control and Computing, 2001. 
[17] M. E. Crovella and A. Bestavros, "Self-similarity in world wide web 

traffic: evidence and possible causes," in ACM SIGMETRICS, 1996. 



Proof: The following set of equations are trivial and 
follow from d23l . 



x (l-r)pf-Ld S = fegL. 

J v{s) J p (j) 



PoCO 

l-poC* 



■ds 



x (l-r)j3T(y) =ln(p (*))|2-ln(l-p (j))|> 
xo(l-r)pT(y)=\n 



Po(0)l-p (y) 



1 



*60 



x (l-r)p 
1 

*o(l-r)j3 



Po(y) i-Po(Q) 
I-P06O Po(0) 

Po(y) 



i - poM 



(38) 



From fl38l ), Po(y) can be expressed as (t24b . ■ 

Proof of Lemma \3\ 

Proof: From ( flOb or (0), practical values of ho must 
satisfy Iiq < 0.5. For ho = and po = 0, the proof is trivial. 
We only need to prove the bounds of jz^ for ho 6 (0,0.5). 

We need (O and the well known property of natural 
logarithms shown next. 



ln.Y < x — 1 



for x > 



(39) 



The proof is by contradiction and has two parts, 
(i) First, we establish that (1 — po)/(l — ho) > 1. Let us 
assume: 



1 -Po 



< 1 



I -h 

From d40l and (fJJ, we get 



ho < Po- 



h < 



In 

- E5 i 

.1-2*0 J 



l-/'o 1 
l-2/7 J 



In 



l-/'o i 
1— 2fto J 



(40) 



(41) 
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Since < h < 0.5, (1 - h Q )/(l - 2h ) > 1 and ln[(l - For arrival rate x E [C/2,°o), since p (0) ^ p* by dHJ) and 
/j ( ))/(l -2/i )] > 0. Multiplying both sides of (Hi} by the p Q (b) := Mo < 26.9% ^ Pg by the Limit property (see Sec. [D, 
denominator term found on the r.h.s. and simplifying, we the critical point exists somewhere y* £ [0,fe]. 



obtain: 



<ln 



I -ho 



1 - 2/i 1 - 2h Q 



(42) 



On the other hand, since \_2h > 0' we can a PPly 
property 



1-2*0 

on . to obtain: 

1-2*0 

, 1 - A 1 - /jo , 
In < 1 = 



Proof of Lemma \5\ 

Proof: From d37] i, it is trivial to see that v(y) is decreasing 
with y since v (y) < 0. Differentiating d37l i and using d45l l. 



l-2h 



o 



1-2/j, 



o 



l-2*o 



(43) 



v"(y)^(l-r) 2 j3 2 p (y) 



which is a contradiction to d42l . Hence, we prove that 
> 1. = 1 when both Mo, &o = 0. 
(ii) Here we establish that (1 — Mo)/(l — *o) < 2. As before, 
let us contradict by assuming that, 



l - Po(y) 
vCv) 



(50) 



1 - Mo 
l-*o 



>2 



Mo < 2Ao-l. 



(44) 



Since v"(y) > 0, v(y) is convex decreasing throughout the 
queue. See Fig. [9b] ■ 

Proof of Lemma [6] 



Proof: The queueing delay is a strict convex increasing 
Since h < 0.5, (04) says that Mo < 0, which is not f unct ion, since from ©5) and ([37), respectively, 



possible. This means the assumption in d44l) must be 
wrong, or that (1 - Mo)/(l - ho) < 2. 
Combining © and (In} completes the proof. 



T (y) = l/v(y) >0 
1 



t 00 



v2(y) 



v(y) 



x (l-r)Ppo(y) 
V-2(y) 



>0 



(51) 
(52) 



Proof of Lemma [5] 

Proo/: Solving for p (y) from (f36t and d37b 

~po(y)-Po(y) 



Po(y)=M 1 ~ r )P 



v(y) 



Proof of Lemma \8\ 

Proof: Set Ar 
(45) Lemma 



Since /3 < 0, and probabilities po(y) < 1, Po(>') is decreasing 
with y. Taking the differentiation further and using d37]>.(l45l 
in place of v (y) and p (y) and simplifying, we obtain 



W)(t(*))|o=i 



t(£>) and a = 1, i.e., ;to2 = *o in 
1 1 



1 +ae^ x o T W/ 3 1 +ae-*oT(*)Mi-i/*) 
1 



Po(y) =jco(i -^)^ 



p (>') v (>')(i- 2 Po()'))-v Cy)Po(y)(i-PoCy)) 



l+fl(l-l//7)-- r T W 

Po(0) 



v2(y) 



2«2po() ; )(i-po()'))(i-3poCy)) 



v2(y) 



(46) 



The critical point (y*,po(y*)) where p (y) = is given as 

1 
3 

1 l-3po(0) 



po(0) + (l-po(0))(l-l/A)"W) 
Po(0)v(0)(l-l/^" T W 

p (0)v(0)(l-l/fe)*o T ( 6 )+ATpi(0)v(0) 
p (b)v(b) 



(53a) 
(53b) 



Po =PoO* 



(47) 



Mo 



(53c) 
(53d) 



1 



v ¥ = - In 
y K 



Po(b)v{b)+N Pl (0)v{0) 
In the above step, we use po(b) = Mo> v(b) = C from ( fTSI l. and 



*2(l-po(0))" 

d48l > is obtained upon substituting p^ for po(y) in d26l >. It 
is easy to see that po(y) decreases in concave fashion on y £ 
[0,y*] and in convex fashion on y £ [y*,£>]. See Fig. [9a] for an 
example. 

Note that the critical point (y*,Po*) exists when the UDP 
arrival rate xq exceeds a certain value x%, calculated using ( fT6b 
as, 



(48) Npi (0)v(0) = (1 - Hq)C from ® . 



Proof of Theorem [2] 

Similar to proof of Lemma [8] we proceed for the general 
a as, 

tt>(T(*)) = ; . t; A^-^wm (54) 



xq = min{x | po(0) ^ Po*} = min 



C I -Mo 
2 1-V 



Above, we use Lemma [3] to state that ^ 1. From 

J, it follows that when xo ^ C/2, po(y) is strictly convex 
decreasing. 



l+a(l-l/b)- ax o< h ) 

But from proof of Lemma [8] above, we find for a = 1 that 
C Mo = 1/(1 - l/b)~ XQZ ^). After rearranging, we obtain 

2- (49) (l-l/fo)-« T w = (l-Mo)/ajUo. Substituting this into ([54]) for 
general a completes the proof. ■ 



